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Abstract 

It is shown that the low-energy single-particle excitation-spectrum of the 
three-band Hubbard model at hole-dopings away from half-filling agrees re- 
markably well with Quantum Monte Carlo data and spectroscopic experi- 
ments within the framework of a conserving approximation that includes self- 
consistently the interaction with charge, spin, and two-particle fluctuations. 
The dispersion of the quasiparticle-peak obtained from the spectral- weight 
function is flat around the (vr, 0) and (0, vr) points as has been observed in re- 
cent angle-resolved photoemission measurement. The significant reduction of 
the quasiparticle-dispersion near the Fermi-energy is due to spin-fluctuations, 

rather than being induced by band effects. 
PACS numbers: 74.20.Mn, 75.10.Lp, Tl.lO.+x 
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The nature of quasiparticle excitations in the high-Tc superconductors (HTSC) is of 
central interest in order to characterize the normal-state of these materials. Angle-resolved 
photoemission spectroscopy (ARPES) shows on the one hand, the existence of a Fermi- 
surface in the normal-state in accordance with Luttinger's theorem supporting the conven- 
tional Fermi- liquid concept. On the other hand, these experiments find a very flat dispersion 
of the quasiparticle excitations near the X=(7r,0) and Y=(0,7r) points. This extended saddle- 
point behavior cannot be derived from single-particle calculations (Hartree-Fock, LDA [^) 
which find, in contrast, a van Hove bandstructure singularity. The importance of many-body 
correlation effects for such a flat quasiparticle dispersion was revealed in recent Quantum 
Monte Carlo (QMC) simulations of both the three-band and the one-band Hubbard 
model [0. Related interpretations were also advanced recently on the basis of QMC sim- 
ulations of the t-J model [||. However, in spite of the success of the numerical techniques 
in dealing with strongly correlated systems, the physical interpretation of purely numerical 
QMC results is often difficult without reference to more transparent theoretical approaches. 

In this letter one-particle spectra for realistic parameters |^,y,|T^] of the three-band Hub- 
bard model are presented based on a conserving approximation |TT]JT^. The obtained spec- 
tral distribution of the quasiparticle peak agrees remarkably well with QMC- and ARPES 
results. The main features of the experimental bandstructure, such as the flat dispersion 
near the Fermi-energy, bandwidth, and Fermi-velocities are reproduced by our theoretical 
approach. Furthermore, in contrast to the QMC simulations, we are able to distinguish the 
influence of different many-body scattering channels. We can clearly show that magnetic 
fluctuations are responsible for the drastic reduction of characteristic energy-scales in the 
quasiparticle excitation spectrum compared to their mean-field values. Thus, the flatness 
has a many-body origin, rather than being induced by band effects, a result which may 
have important implications particularly for mechanisms of superconductivity based on spin 
fluctuations |r3|-p^. The results presented here were obtained at much lower temperatures 



than in numerical simulations such that a direct comparison with experiments is possible. 
Our analysis is based on the two-dimensional three-band Hubbard model 0,01 with 



the Hamiltonian H — = Hq — fiN + Hu, where 

Hq- fiN = ^(^d - fi)ndia + ^{£p- f^)npja 
+ H Uj{dl^pja + h.c.) , and 

<i,j>,<j 

Hu = UdJ2 ^di^^dii ■ (1) 
i 

Here and pj^. denote the electron creation operators for Cu-3dx2_y2 and 0-2p orbitals at 
sites i and j, respectively with spin a, and ridia = d\^di„. The one-electron part Hq — fiN 
is determined by the local orbital levels Sd (Cu) and Ep (O) with charge-transfer energy 
A = Ep — Ed, and a Cu-0 hopping tij = ±t between nearest-neighbor Cu-0 sites < i,j >, 
where the sign takes into account the phase factors for the 2p and 3d^.2_j^2 orbitals. fi is the 
chemical potential and the interaction part Hu consists of the on-site Coulomb repulsion 
Ud (Cu). The calculations were performed for the parameter set Ud = Q, A = 4 in the hole 
representation (as everywhere, the energy-unit is t = 1). This choice follows from recent 
QMC studies ||^J^,|1^ , where it was shown that the above parameter set leads to a consistent 
description of several important features of the HTSC. 

Bickers and Scalapino have developed a propagator-renormalized fluctuation-exchange 
(FLEX) approximation |]TT| for the one-band Hubbard model based on the conserving ap- 
proximation (CA) of Baym and Kadanoff []TB[ which has resulted in significantly improved 
agreement with QMC results |T^. Properties of the superconducting state in the same model 
were also studied recently by that method The starting-point of this self-consistent field 



solution is the non-interacting Hamiltonian Hq in Eq. |l|. Diagonalization of Hq leads to the 
eigenvalues i?°(k) (n=l,2,3) which build up the unperturbed bandstructure. The single- 
particle propagator of the non-interacting system can now be expressed in terms of these 
energies £'°(k) and the corresponding eigenvectors c„i,(k): 

"^0/1^ , , \] — sr^ c„^(k)c„^/(k) 

which becomes a matrix labeled by the orbital indices (z/= d, Px, Py) of the three-band 
Hubbard model. We employ the Matsubara (T ^ 0) formalism where the discrete fermionic 



frequencies are given by cUm — (2m + l)7rT , m — {..., —2, —1, 0, 1, 2, ...}. The Coulomb 
interaction Hu is accounted for within the field-theoretical FLEX approximation by a power 
series for the self-energy This series includes diagrams (symmetric in particle-hole and 
particle-particle channels) taking into account the interaction of electrons with density (d), 
spin (s), and particle-particle (p) fluctuations. For the present three-orbital Hubbard model 
the single-particle self-energy is as follows: 

nk)\^,= '^T.{ [G{k-q)\^^{V^'\q) + V^'\q) + V^^\q)) + [G{-k + q)\^V^^\q) } , 

(3) 

where 



v^'\q)^ul-xA<i)^ 

V^^\q) = -Ulxpp{q) 



1 



- 1 



(4) 



,1 + UdXppiq) 

k stands for the fermionic momenta and Matsubara frequencies (k, cUm), and q for the bosonic 
(q, Un) ones. Xph and Xpp are the propagators of non-interacting particle-hole and particle- 
particle fluctuations: 

-T 



Xph{q)-^i:[G{k + q)]^^ [G{k) 



Xpp{q)-i;7i:[Gik + q)]^^[G{-k) 



dd ' 



dd 



(5) 



Inserting this self-energy into Dyson's equation the single-particle propagator has to be 
computed self-consistently: 



G'(k, UJm) 



G°(k,a;jl - \j:[G,U](k,u;m) 



(6) 



with G° standing for the single-particle propagator of the non-interacting system. 
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The single-particle excitations are determined by the peaks of the spectral- weight func- 
tion 

^(k,o;) = --^ ^[^(k,^;)]^^ . (7) 

For the analytic continuation of the Green's function G to real frequencies lo we employ a 
simple Fade algorithm |2^. All results were obtained on k-space lattices of 16 x 16 points 
in the first Brillouin zone (BZ). 

In Fig. [l|, we plot the dispersion of the quasiparticle peak in A(k, cu), which is located 
closest to the chemical potential yU along the symmetry lines in the BZ. The hole doping 
[b = 0.25 away from half-filling) and the relatively high temperature (T = O.lt) were chosen 
to compare our FLEX approximation to our best converged QMC simulations which were 
performed on a 8 x 8 lattice (192 sites) with periodic boundary conditions. The FLEX 
dispersion for the single-particle excitations is in good agreement with the QMC data, in 
particular, for the peaks in y4(k, cj) crossing the chemical potential /i (cj = in Fig. |1|). 

In order to test the influence of the temperature on the spectrum, we have also performed 
both FLEX calculations and QMC simulations for the even higher T = 0.5t for various 
dopings. Like in similar investigations for the single-band Hubbard model at the electronic 
filling (n) = 0.83 at the same temperature [^], we find the u = crossing in A(k, u) slightly 
above the fiat band region at X = (vr,0), whereas this crossing occurs below this region at 
T = O.lt. This indicates a strong dependence of the a; = crossing on T at least in the 
high-temperature regime. 

For the much lower temperature T = O.Olt ~ 150K, Fig. |^ displays the spectral weight 
A{]s.,uj) for energies near the chemical potential /i and the three k-space lines as in Fig. 
|1]. The hole doping was set to 6 = 0.15 away from half-filling in order to compare with 
spectroscopic measurements in the metallic phase of the HTSC. The pairs of numbers on 
the right side of each panel have to be multiplied by 7r/8 to give the corresponding k- vectors 
in the two-dimensional BZ. These curves can be directly compared with experimental data 
obtained from angle-resolved photoemission spectroscopy (ARPES) PJ^. Good agreement 



is found for the quasiparticle excitation spectra calculated by the FLEX approximation and 
the spectroscopic measurements. The main features of the experimental data are reproduced 
by our theoretical approach: Along the F — M high symmetry direction we find a dispersing 
quasiparticle peak which crosses through the Fermi-energy near k = (7r/2,7r/2). The two 
other symmetry lines in Fig. ^ clearly show the existence of the flat band dispersion near 
X=(7r,0), the most striking result of the ARPES experiments. The quasiparticle peak of 
A{]s.,uj) is broadened for excitation energies away from fi due to strong correlation effects. 
In fact, it can be seen that not only the position of the maxima in A{]<.,uj) (dispersion) but 
also the overall shape of the spectra agrees with the experimental results. 

The single-particle bandstructure obtained from the low- lying peak in A(k, u) is depicted 
in Fig. 1^ in comparison to results of ARPES experiments in Bi2212. In order to discern the 
influence of different many-body scattering channels, we present the data of three approxi- 
mative stages: In the Hartree-Fock (HF) or band-theory approximation (dashed line in Fig. 
^ the HF bandwidth, Fermi- velocities, and the characteristic energy-dispersion at (7r,0) are 
far away from agreement with the experiment. Charge fiuctations given by \/W(g) in Eqs. 
and ^ tend to reduce these differences (dotted line in Fig. ^, but only the inclusion of 
both longitudinal and transverse magnetic fluctuations |]22| V^^\q) leads to a flat quasi- 
particle band near the Fermi-energy. The particle-particle channel V^^\q) contributes only 
in a minor way to the full FLEX calculations (full line in Fig. |]) for the given parameter 
set. Therefore, Fig. ^ demonstrates that spin- fluctuations are responsible for the reduced 
energy-scales compared to mean-field calculations, especially for the flat dispersion near n 
and the reduced saddle-point energy at (7r,0). The most dominant contributions stem from 
magnetic modes with a wave- vector q near to the antiferromagnetic point (vr,7r). 

Finally, Fig. 4 shows a strong correlation between the energy difference A = E{X) — n 
and the energy cUgpin of the maximum in ImXspin(q, cu) at the incommensurate wavevector, 
where Imxspin reaches its largest value in the BZ. This suggests that for dopings 5^0.25, 



the energy difference A is locked-in to the magnetic excitations close to the AF-point [C 



Moreover, we observe in our FLEX-calculations a linear tu-dependence of Im S down to 



approximately the same energy-scale. Below this energy-scale, Im S ~ cu^ as expected for a 



Fermi-liquid. Details of these features will be given elsewhere p^j. 

In summary, we presented the single-particle excitation spectrum of the three-band Hub- 
bard model within the framework of a conserving approximation. This many-body approach 
is based on a self-consistent evaluation of Dyson's equation including fluctuation-exchange 
(particle-hole and particle-particle) self-energy contributions. The results are in remarkable 
agreement with Quantum Monte Carlo data and spectroscopic experiments for low-lying 
single-particle excitation energies. The low-energy physics of this FLEX calculation is con- 
trolled by many-body fluctuations in the magnetic scattering channel leading to a significant 
reduction of characteristic bandstructure energy-scales from their mean-field values. 
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FIGURES 

FIG. 1. Quasiparticle dispersion along three cuts through the Brillouin zone with the 

high-symmetry points r=(0,0), X=(7r,0), M=(7r, vr). The FLEX result for hole-doping (5=0.25 
and temperature T = 0,lt is compared with QMC data. 

FIG. 2. Spectral- weight function A{li,Lo) for excitation energies near the chemical potential 
and k- vectors along high-symmetry lines in the BZ. 

FIG. 3. Quasiparticle dispersion along high-symmetry lines in the BZ. Hartree-Fock (HF) re- 
sults (dashed line) and FLEX calculations which include only the charge scattering channel (dotted 
line) and all symmetrically chosen particle-hole and particle-particle channels (full line) are com- 
pared to ARPES experiments (open circles: Ref. [2], full circles: Ref. [3]). 

FIG. 4. Energy-difference E{X) — fi and energy of magnetic excitations ct'spin 

close to the 

AF-point as a function of doping. 



10 



r — >x 



X — > M 



r — > M 




-0.4 0.0 0.4 -0.4 0.0 0.4 -0.4 0.0 0.4 



Energy relative to the chemical potential ji 



1.0 
0.75 

0.5 
0.25 

0,0 
-0.25 
-0.5 
-0.75 
-1.0 
-1.25 



_ Ud = 6 




N 

S 


A = 4 


y 

/ 

/ 


\ 

\ 

\ 


~ 6 = 0.15 


/ 

/ 

/ 

/ 

/ 


\ 

\ 

\ 

•.. \ 


- /? = 100 


/ - — 

/ 

/ 


\ 

^\ • \ 
\ . \ 

N. ■ • ^ 
-A 
















all fl. 






\ ■•. 
\ ■ •• 
\ ■• 


/ 


charge fl. 


\ 

\ 

\ 




\ 




HF 


\ 

\ 

\ 


/ 


• Exp.l 


\ 

\ 

\ 

\ 


/ 

/ 


o Exp. 2 


\ 

\ 

\ 






\ 



r 



X 



M 



r 



